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Carrier-induced ferromagnetism is investigated in a diluted Hubbard model for ordered impurity 
arrangements. The delicate competition between particle-hole processes contributing to the spin 
couplings results in a rich variety of behaviour. The ferromagnetic transition temperature obtained 
within the spin-fluctuation theory is in good agreement with reported values for Gai_ x Mn x N. 



I. INTRODUCTION 

The recent discovery of ferromagnetism in diluted mag- 
netic semiconductors (DMS) such as Gai_ x Mn x Asi*^ 
and Ga!_ x Mn x N,- and the intensive efforts to in- 
crease the ferromagnetic transition temperature in view 
of potential technological applications, has generated 
tremendous interest in the novel ferromagnetism ex- 
hibited by these systems in which magnetic interac- 
tion between localized spins is mediated by doped 
carriersi 4 ! 5 ! 6 !!! 8 ! 9 ! 1 "! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 1 ?, 1 ^^ The nature of 
the ferromagnetic state has also attracted much atten- 
tion, particularly in the context of competing antiferro- 
magnetic (AF) interaction which result in noncollinear 
ordering, spin-glass behaviour, and significant sensitivity 
of spin stiffness and transition temperature T c on carrier 
concentration, spin clustering etc. 

DMS are mixed spin-fermion systems, involving ran- 
domly distributed localized magnetic impurities (e.g., 
Mn ++ , having S — 5/2) and mobile carriers (e.g., holes) 
in the semiconductor band. With carrier concentration 
p much smaller than the magnetic impurity concentra- 
tion x, the DMS systems provide a complimentary limit 
to Kondo systems. Conventionally, the coupling between 
localized impurity spin S and mobile valence band holes 
is represented by the exchange interaction — JS.cr, where 
er = ^[cr] 1 ]/ is the fermion spin operator. 

Cuprates form another class of strongly correlated sys- 
tems in which the concentration of doped carriers is small 
in comparison to that of (Cu) spins. Extensively studied 
within the Hubbard model, the dominant interaction in 
cuprates is the AF exchange interaction between neigh- 
bouring Cu spins, and AF spin correlations persist even 
for small hole doping. Against this strong tendency for 
AF ordering, the delocalization energy gain of doped car- 
riers, which favours parallel spin alignment, results in 
only a marginal twisting of spins, as in the proposed spi- 
ral phasesiS2i2i It is therefore interesting that elimination 
of the strong AF spin interaction in a diluted Hubbard 
model^i with no direct hopping between relatively dis- 
tant impurity spins, does indeed lead to a ferromagnetic 
state stabilized by carrier-induced spin couplings. 

The spin stiffness in the ferromagnetic state of the di- 
luted Hubbard model goes through a maximum with re- 
spect to both doping concentration p and the interaction 
strength U. This optimization behaviour can be qualita- 
tively understood in terms of a competition between the 
increasing magnitude of carrier spin polarization x{U) 



and the increasing rapidity of its oscillation, which limits 
the growth of the spin coupling U 2 Xij(U) between two 
magnetic impurities at a fixed separation. Similar be- 
haviour was observed in the spin-fermion model for the 
effective ferromagnetic coupling = J 2 Xij(J) m terms 
of the generalized spin response x(J) f° r finite JH As 
X(J) involves particle-hole processes, the spin couplings 
identically vanish in the absence of doping. However, in 
the diluted Hubbard model, even with no carrier doping 
(p = 0) in the majority-spin band, particle-hole processes 
involving the empty minority-spin impurity band result 
in antiferromagnetic spin couplings, which destabilize the 
ferromagnetic state below a critical doping concentration. 

In this paper, we study the novel ferromagnetic state 
of the diluted Hubbard model for ordered impurity ar- 
rangements. We focus on magnon excitations and the 
spin stiffness, which provide quantitative measures of the 
stability and spin couplings in the ferromagnetic state, as 
well as the transition temperature in three dimensions. 
While preliminary real-space studies for both ordered and 
disordered impurity arrangements were carried out ear- 
lier on three-dimensional finite-size lattices^ use of k- 
spacc representation in this paper allows for much larger 
lattices, thus permitting a more refined study of the com- 
peting spin interactions with respect to carrier concen- 
tration, impurity separation, interaction strength, and 
wave vector. Indeed, the antiferromagnetic - ferromag- 
netic quantum phase transition stands out as a signifi- 
cantly prominent feature in the low doping regime. AF- 
F and spin-glass transitions at low doping have actually 
been observed in Mn-doped II- VI semiconductors. 22 - 23 - 24 
With increasing doping concentration, an anomalous in- 
crease in the spin stiffness is observed which involves, as 
explained in section III B, a subtle interplay of impurity 
moment reduction, impurity character of doped states, 
and competing particle-hole processes, all characteristics 
of the itinerant ferromagnetic state. 

Within the relatively simpler ferromagnetic Kondo lat- 
tice model (FKLM), in which a localized spin is present 
at every lattice site, magnon excitations have been stud- 
ied recently, as a function of electron density n in the 
conduction band and the spin-fermion coupling J, in the 
context of heavy fermion materials^ ferromagnetic met- 
als Gd, Tb, Dy, doped EuX, 26 and manganitesi 27 i 28 i 29 
Magnon dispersion has also been obtained in the con- 
text of DMS i 30 i 31 ' 32 i 33 However, a uniform impurity- 
induced Zeeman splitting of the carrier spin bands is 
assumedi 30 ' 31 ! 33 
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FIG. 1: An ordered arrangement of impurity atoms (•) on a 
square host lattice (o). Also shown are the sublattice labels 
corresponding to the four distinct sites within the unit cell. 



II. DILUTED HUBBARD MODEL 

We consider a diluted Hubbard model 

H = t~^2 a i t a a i+S,a +t' ^ a\ a ai+8,v 

+ e d a \c a i<y + u ( n n - ( n i)) ( n n - ( n i)) 

1,0 I 

(1) 

on square and cubic lattices with nearest-neighbour (NN) 
hopping between sites i and i + 5. Here / refers to the 
impurity sites, ed is the impurity on-site energy, and 
(ni) = (nii + n il)/2 is the spin-averaged impurity den- 
sity. The energy-scale origin is set so that the host on-site 
energy is zero, and we take the impurity level to lie at the 
top of the host band (e^ = 2dt = W/2 in d dimensions) to 
optimize local-moment formation. For a positive sign of 
the hopping term, the top of the host band lies at k = 0, 
as in the valence band of DMS. For simplicity, we take the 
same hopping (t' = t) between the host-host and host- 
impurity pairs of sites. Higher spin magnetic impurities, 
such as the S — 5/2 Mn impurities in Gai_ x Mn x As, can 
be realistically represented within a generalized Hubbard 
model involving multiple orbitals per site^i 



A. Hartree-Fock ferromagnetic state 

In the Hartree-Fock (HF) approximation, the interac- 
tion term reduces to a magnetic coupling 



H; 



HF 



I 



a-i.A] 



(2) 



with the local mean field A/, resembling the semi- 
classical limit of the exchange interaction — Y] t Jer/.S/ 
in the spin-fermion model. Here the impurity spin opera- 
tor (j i = fyj[(r]^>i in terms of the fermionic field operator 

9 j = ^ ^ . The mean field A/ is self-consistently de- 
termined from the ground-state expectation value: 

2A, = (crj)U . (3) 
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FIG. 2: Variation of impurity magnetization with hole con- 
centration for the two-dimensional case (x = 25%). With 
increasing U, the impurity character of majority-spin states 
at the top of the band diminishes, and the impurity moment 
is therefore less susceptible to hole doping. 



B. Sublattice-basis representation 

We consider ordered (superlattice) arrangements of 
magnetic impurities on square and cubic host lattices, 
with impurity separations of 2a and 3a. Translational 
symmetry within the sublattice basis conveniently allows 
Fourier transformation to momentum space. For con- 
creteness, we consider a square host lattice in the follow- 
ing, with magnetic impurities placed at every other host 
site, corresponding to superlattice spacing 2a and impu- 
rity concentration x = 25%. There are four sublattices, 
numbered a = 1,2,3,4, corresponding to the four sites 
in the unit cell, as shown in Fig. 1. We choose length 
and energy units such that the lattice spacing a = 1 and 
the hopping term t = 1. 

Without loss of generality, we assume a mean field 
A/ = Az in the z direction. Fourier transformation 
within the four-sublattice basis yields the HF Hamilto- 
nian 
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for spin cr, where e£ = 2tcosk x and = 2tcosk y corre- 
spond to hopping terms in the x and y directions, respec- 
tively. Here the field operator 'Fka = (a\. a a\ a a\ a a^ a ) 
defines the sublattice basis, where ag CT refers to the 
fermion operator for sublattice index a. Generaliza- 
tion to other impurity concentrations and dimensions is 
straightforward. Hopping terms to the right (up) and left 
(down) do not connect the same sublattice in general, 
yielding a Hcrmitian Hamiltonian matrix with hopping 
terms te ±lkm etc. The case of impurity spacing 2 is, how- 
ever, special as hopping terms in opposite directions do 
connect the same sublattice, yielding a real symmetric 
Hamiltonian matrix. 

The HF Hamiltonian matrix is numerically diagonal- 
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FIG. 3: Total and impurity (I) density of states for the 2-d 
case (x — 25%) for two values of the mean field A, showing the 
split-off impurity band and the impurity-induced sub-bands. 

ized to obtain the four eigenvalues E£ a , corresponding 
to the four sub-bands fi — 1,2,3,4 in increasing order 
of energy. The four-component eigenvectors <j>^ yield 
the amplitude on sublattice a. Summing over occupied 
states yields the impurity magnetization and the self- 
consistency condition 

^=m(A) = (af)= £ (C =1 ) 2 -(C =1 ) 2 - ( 5 ) 

Variation of the impurity magnetization m with doping 
concentration p is shown in Fig. 2. 

C. Quasiparticle spectrum 

Each magnetic impurity contributes two spin states. 
The minority-spin impurity state (energy a + A) is 
pushed up by the local Coulomb repulsion, and forms 
a split-off narrow impurity band due to the small overlap 
with neighbouring impurity states. The majority-spin 
impurity state (energy e<j — A) mixes with the host band 
states for — A <C W, and splits off on the low-energy 
side for large A. Figures 3 and 4 show a comparison of 
the quasiparticle density of states N(E) for two A val- 
ues for the 2-d and 3-d cases, showing the four and eight 
sub-bands corresponding to the sublattices, respectively. 



FIG. 4: Density of states for the 3-d case (x — 12.5%) for two 
values of the mean held A. 

For minority spin, while the impurity band shifts up 
with increasing A and narrows due to decreasing effective 
impurity hopping t c g ~ t 2 /(e^ + A), the host sub-bands 
remain unaffected reflecting strong host-impurity decou- 
pling. The majority-spin sub-bands, however, are pulled 
down in energy by the exchange interaction, decreasing 
the energy separation between the two bands. With in- 
creasing A, the impurity character of majority-spin states 
at the top of the band diminishes, as seen in Fig. 3. 

In the undoped insulating state, the Fermi energy lies 
in the gap between the top of the majority-spin valence 
band and the minority-spin impurity band. The empty 
impurity band ensures local-moment formation on each 
impurity site. Carrier doping is introduced by adding 
holes to the top of the majority-spin band. The mean 
field A decreases with doping concentration p, which 
shifts the majority-spin band to the right. The Fermi 
energy therefore remains above the minority-spin valence 
sub-band (/z =3 or 7), yielding a hnite Stoner gap for low 
doping. 

III. TRANSVERSE SPIN FLUCTUATIONS 

Transverse spin fluctuations are gapless, low-energy 
excitations in the broken-symmetry state of magnetic 
systems possessing continuous spin-rotation symmetry. 
Therefore, at low temperatures they play an important 
role in diverse macroscopic properties such as existence 



of long-range order, spin stiffness magnitude and temper- 
ature dependence of magnetization, transition tempera- 
ture, spin correlation etc. 

We study the time-ordered, transverse spin propagator 



xt j -(t-t') = i(* G \T[S+(t)Sr(t>m G ) 



(6) 



in the the ferromagnetic ground state \^g), involving 
the spin-raising (S^~) and spin-lowering {Sj) operators 
at sites i and j. In the random phase approximation 
(RPA), the magnon propagator 



[x + -(q^)] = [x- + (-q,-")] = r 



[X°(q,")] 



W(q,")] 



(7) 



in q, lu space, where the zeroth-order, particle-hole prop- 
agator [x°(q, u))] is obtained by integrating out the 
fermions in the self-consistent broken-symmetry state. In 
the sublattice basis we have 

[X°(q,w)]a/J = i f ^^[&(Ku')U[G l (y,u>' + u)]{, 
J n k 

E k'i >Ep J,M« fh^f*^! 3 fh^a* 



E»<E F k ^ 



E^>E F 

£ 

E» kll <E F 



^kj 



(8) 



where k' = k + q. In Eq. (7), the diagonal interaction 
matrix [U] aa — US al has non-zero element only on the 
impurity sublattice a = 1. It is therefore convenient to 
write 



[x + -(q, 
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1 



(9) 



in terms of a matrix L4(q,w)] = [U] - [U]{x°(q_, w)][U], 
which has non- vanishing matrix elements only on the im- 
purity sublattice. Magnon-mode energies w q are there- 
fore obtained from the magnon pole condition 



l-tfX°(q,w q ) 







(10) 



where \° (q> w ) represents the impurity-sublattice matrix 
element of [x°(q> w )]a/3- I n the ferromagnetic state, typ- 
ically x°(q, w) = jj — Aq 2 + Blo for small q,ui, so that 
the magnon energy w q = (A/B)q 2 . The coefficient A is a 
measure of the spin stiffness, whose sign determines the 
stability of the HF ferromagnetic state with respect to 
transverse fluctuations, as discussed below. 



A. Stability of the ferromagnetic state 

The undoped state, with filled majority-spin states and 
an empty minority-spin impurity band, amounts to a 
half-filled case with respect to magnetic impurities. Vir- 
tual hopping between impurity sites therefore generates 
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FIG. 5: (a) Negative magnon energy in the absence of hole 
doping indicates instability of the ferromagnetic state due to 
antiferromagnetic spin couplings, which weaken with increas- 
ing U. (b) Stabilization of the ferromagnetic state with in- 
creasing hole concentration. 



antiferromagnetic exchange interaction t 2 s /U between 
neighbouring magnetic impurities, resulting in antiferro- 
magnetic ordering. The carrier-mediated ferromagnetic 
couplings are absent, and the self-consistent ferromag- 
netic state is therefore unstable and actually represents 
a saddle point. The instability is reflected in a negative 
sign of the coefficient A (spin stiffness) , which makes the 
transverse response eigenvalue E/x°(q) > 1 and yields 
negative magnon energies [Fig. 5(a)]. With hole dop- 
ing, ferromagnetic coupling strengthens, long- wavelength 
magnon-mode energy changes sign, and ferromagnetic or- 
dering is stabilized at some critical hole concentration 
[Fig. 5(b)]. Thus the diluted Hubbard model is charac- 
terized by an AF - F transition at a finite doping concen- 
tration, which is more prominent for small U. 

A significant feature of Fig. 5(b) is the distinct be- 
haviour of long- and short-wavelength magnon modes. 
Even in the unstable regime, with negative-energy long- 
wavelength modes, the zone-edge modes (q x ,q v ~ t/2) 
may have positive energy, thereby giving a spurious in- 
dication of stability. Incorporating only the zone-edge 
(Ising) excitations, the dynamical mean-field theory is 
therefore insensitive to the non-trivial long-wavelength 
behaviour arising from competing spin couplings. 
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FIG. 6: Variation of spin stiffness with (a) hole concentration FIG. 7: Variation of spin stiffness with (a) hole concentration 
and (b) interaction strength U, for the 2-d case (x = 25%). and (b) interaction strength U, for the 3-d case (x = 12.5%). 



B. Spin stiffness 

The long-wavelength magnon-mode energy w q and 
spin stiffness uj^/q 2 provide a composite measure of impu- 
rity spin couplings in the carrier-mediated ferromagnetic 
state. Figures 6 and 7 show the behaviour of spin stiffness 
with doping concentration p and interaction strength U 
for the 2-d and 3-d cases with impurity spacing 2. A 
minimum hole concentration to stabilize the ferromag- 
netic state is clearly seen, especially for low U, when the 
competing antiferromagnctic spin couplings are relatively 
stronger. The increase in spin stiffness with hole con- 
centration exhibits the essence of carrier-mediated ferro- 
magnetism. The initial increase is nearly linear in two 
dimensions and distinctly sublinear in three dimensions. 
The spin stiffness exhibits an optimization with respect 
to both p and U , which can be understood qualitatively 
within the generalized RKKY theory in terms of a com- 
petition between increasing magnitude of carrier-spin po- 
larization and increasing rapidity of its oscillation. 

Similar behaviour of spin stiffness is seen (Fig. 8) when 
impurity spacing is increased to 3 in the 2-d case. The 
peak spin stiffness is clearly reduced, reflecting the weak- 
ening of spin couplings with increased impurity separa- 
tion. However, in the low doping regime, the spin stiff- 
ness for x = 1/9 is actually greater than that for x = 1/4 
(see inset). As the stiffness must eventually decrease with 
increasing impurity separation, this interestingly shows 
an optimization behaviour with respect to impurity con- 



centration x as well. 

The anomalous increase at higher p, seen in Figs. 6- 

8 for lower U values, is due to a competition between 
two particle-hole processes, as discussed below. Figure 

9 shows the doping dependence of U 2 [x°(0) — x°( c l)]/9 2 
for the two particle-hole processes in Eq. (8), labeled as 
1 and 2, which provide measures of the corresponding 
spin-stiffness contributions. Positive (negative) value in- 
dicates dominant ferromagnetic (antiferromagnctic) spin 
coupling. Process 2, which is activated only upon dop- 
ing in the majority-spin band, is seen to be dominantly 
ferromagnetic in the low-doping regime and exhibits the 
typical RKKY peak behaviour before turning negative 
due to the oscillating carrier-spin polarization. However, 
process 1, which is weakly antiferromagnetic for p w 0, 
changes sign with doping, and becomes increasingly fer- 
romagnetic; it is the sharp increase seen at higher doping 
which causes the distinct minimum feature in the total 
spin stiffness. This feature is more pronounced at lower 
U because the sharply reduced impurity moment and A 
at higher doping (see Fig. 2), cause the impurity states to 
be more extended, thus enhancing the process 1 stiffness. 



C. Transition temperature 

Determination of the magnon spectrum in the carrier- 
mediated ferromagnetic state allows an estimation of the 
transition temperature T c for the diluted Hubbard model 
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FIG. 8: Increasing the impurity spacing to 3 results in qual- 
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FIG. 9: Spin-stiffness contributions corresponding to the two 
particle-hole processes in Eq. (8), showing the distinct min- 
imum feature and the anomalous increase arising from their 
competition. Here U — 9. 



in three dimensions. As the ferromagnetic state is charac- 
terized by small spin stiffness due to strong competition 
between ferromagnetic and antiferromagnetic spin cou- 
plings, the dominant contribution to reduction in magne- 
tization is from the thermal excitation of long- wavelength 
magnon modes. Therefore, the transition temperature 
T c can be estimated in terms of an equivalent Heisenberg 
model with matching spin stiffness. 

For a nearest-neighbour, spin-5 quantum Heisenberg 
model with interaction energy J on a cubic lattice (coor- 
dination number z = 6), the magnon energy 



D q = Z,JS{\ - 7 q ), 



(11) 



where 7 q = (cos 2q x + cos 2q y + cos 2q z )/3, corresponding 
to the magnetic lattice spacing 2. Indeed, the magnon 
energy is maximum for q x = q y = q z = tt/2, as also seen 
in Fig. 5 for the 2-d case. Considering the small q limit 
of Eq. (11), and matching the spin stiffness, we obtain 
D = ujq/q 2 = iJS = 2 J, where we have set 5=1/2 



FIG. 10: Comparison of the transition temperature T c ob- 
tained from the renormalized spin-fluctuation theory (SF) 
with the HF result (temperature scale on the right), high- 
lighting the order-of-magnitude difference between the spin- 
ordering and moment-melting temperatures. Here x = 12.5%. 



for the diluted Hubbard model having a single magnetic 
orbital per site. 

For the spin-5 Heisenberg model, the transition tem- 
perature 



zJS(S+l) 



f£ = DS(S+l)f£ 



(12) 



within the renormalized spin-fluctuation theory is some- 
what lower than the mean-field value T C MF = zJS(S + 
l)/3 due to the spin-fluctuation factor Jsf = S q (l — 
7q)~\ which is approximately 1.5 for the cubic lattice. 
Extrapolating to spin-5 magnetic impurities in the di- 
luted Hubbard model, represented by multiple magnetic 
orbitals per site^ with the same equivalent interaction 
energy J = D/2 as obtained above, the transition tem- 
perature reduces to the expression in Eq. (12). Taking 
5 = 5/2 for Mn impurities, and a realistic heavy-hole 
bandwidth W = I2t = 2 eV, the transition temperature 
evaluated from Eq. (12) is shown in Fig. 8. Also shown 
for comparison is the HF result obtained from a finite- 
temperature analysis of the self-consistency condition. 
The HF T c is naturally an order of magnitude higher as 
it really corresponds to a moment-melting temperature. 

At the HF level, the impurity moment vanishes at a 
temperature T ~ A, the minority-spin impurity-band en- 
ergy relative to the Fermi energy. The marginal increase 
seen in T C HF with doping is simply due to the decreasing 
Fermi energy, whereas the impurity-band energy remains 
essentially unchanged. 

Similar estimations of T c from magnon energy were 
recently carried out for the ferromagnetic Kondo lattice 
model^ and in the context of DMS; 17 i 30 i 31 ' 33 as a func- 
tion of carrier concentration p. Optimization behaviour 
was found in these cases as well, with T c increasing with 
p upto a maximum, followed by a monotonic decrease. 
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IV. CONCLUSIONS 

Carrier-induced ferromagnetism was investigated in 
the diluted Hubbard model for ordered impurity arrange- 
ments. Momentum-space representation within the sub- 
lattice basis permitted a detailed study with respect to 
carrier doping, impurity spacing, electron correlation, 
and wave vector. The delicate competition between spin 
couplings results in an AF-F quantum phase transition in 
the low-doping regime. Competition between the increas- 
ing magnitude and increasing rapidity of oscillation of the 
carrier-spin polarization was observed, yielding a char- 
acteristic optimization of the spin stiffness with doping 
concentration p and interaction strength U . Surprisingly, 
in the low doping regime, the spin stiffness was found to 
actually increase when the impurity separation increased 
from 2 to 3 in the 2-d case, indicating an interesting opti- 
mization with respect to impurity concentration as well. 
In addition, the itinerant ferromagnetic state exhibits a 
subtle interplay of impurity moment reduction, impurity 
character of doped states, and competing particle-hole 
processes in the carrier-spin polarization, resulting in a 
distinct minimum and an anomalous increase in the spin 
stiffness at higher doping concentration. 

For the cubic impurity arrangement, doping depen- 
dence of the ferromagnetic transition temperature T c 
was calculated within the renormalized spin-fluctuation 



theory from the magnon spectrum, and compared with 
the HF result so as to highlight the difference between 
spin-ordering and moment-melting temperatures. For 
5 = 5/2 and a realistic heavy- hole bandwidth of 2 eV, 
the peak calculated T c of 960 K at p sw 3% is in good 
agreement with the highest reported T c value of 940 K 
observed in Gai_ x Mn x Ni£ 

While our main objective was to examine the novel fer- 
romagnetism in the diluted Hubbard model, several fea- 
tures of this itinerant ferromagnetic state are of relevance 
to DMS systems as well. For example, impurity bands 
formed by overlap of Mn d orbitals are prominent fea- 
tures in density functional calculations within the local 
spin density approximation (LSDA) , and a characteristic 
dependence of T c on Mn concentration is obtained for 
impurity bands formed within the gap, as is the case for 
Gai_ x Mn x Ni^ Furthermore, correlation effects have also 
been examined recently within the LSDA+U approach 36 
to obtain better agreement of the Mn d-DOS peak posi- 
tion with recent photoemission experiments. 
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